Astronomy & Astrophysics manuscript no. cbc' low latency pipeline 


©ESO 2012 


January 16, 2012 





First Low-Latency LIGO+Virgo Search for Binary Inspirals and their 

Electromagnetic Counterparts 

J. A badiffl B. P. Abbot0, R. Abbot0, T. D. AbbotPI M. AbernathjH! T. Accadic0, F. AcernesfiJ^ C. AdamJH, R. AdhikarR C. AffeldCH M. AgathoJi, K. AgatsumJH, R Ajitld, 
B. AlleiEUm E. Amador Ceroiin! D. AmariuteP, S. B. AndersoiP, W. G. AndersoiP, K. AraR M. A. AraiiP, M. C. ArayiP, S. M. AstoiP, R AstonJHJ, D. AtkinsoiP, R AufmutlEEl 
C. AulberOU B. E. AylotP, S. BabalP, R BaketHI G. BallardiiP, S. BallmeP, J. C. B. BaragoyJH D. BarkeP, R BaronfiJ^ B. Bar0, L. BarsottP, M. BarsugliP, M. A. BartoiP, 
I. BartoP, R. Bassir0, M. BastarrikiCl A. Bast&, J. BatclP, J. BauchrowitJzH Th. S. BaueP, M. BebronnJi! D. BeclP, B. BehnkJH, M. BejgeilJ, M.G. BekefiJ, A. S. Bel0, 
A. BelletoilJ^ I. BelopolskP, M. Benacquist<P, J. M. Berlinefni A. BertolinCEl J. Betzwiese0, N. Beveridgfl R T. BeyersdorP, I. A. BilenkJU, G. Billingsle^, J. Bircl0, 
R. BiswaJH, M. BitossP, M. A. BizouarcC^}, E. Blacfl J. K. Blackbun0, L. BlackburiP, D. BlaiEU B. BlancCirM. BlonP, O. BoclJzH T. R BodiycP, C. BogaijZiJ R. BondaresciP, 
R BondiP', L. BonellEl^ R. Bonnan(!HI R. Borfl M. BoyM\ V. BoschP, S. BosJH, L. BosEl, B. BouhojSl S. BraccinP, C. Bradaschid, R R. BradjEH, V. B. BraginskjUl 
M. BranchesP^, J. E. BraJH J. BreyeQS T. BrianP, D. O. Bridget, A. BrilleP, M. BrinkmaniQEI V. BrissoiHJ, M. BritzgeiB A. R Brook0, D. A. BrowiP, T. BulilB^, 
H. J. BulteI0^ A. BuonannJH, J. Burguet-CastelP, D. Buskuli(0, C. BujEH, R. L. ByelH L. CadonatEH G. CagnolSl, E. Callon0^ J. B. CamfH^ R Campsifl J. Cannizz(P, 
K. CannoiP, B. CanueP, J. CacP, C. D. Capan(P, R CarbognanP, L. CarboniHl S. CaridJ^, S. CaudilP, M. Cavagli^ll R CavalieJU, R. CavalierP, G. Cell<&, C. Cepedcfl 
,^ E. CesarinP, O. ChaibP, T. Chalermsongsafl R CharltoiE^ E. Chassande-MottiiCH S. ChelkowskP, W. CheiP, X. CheiP, Y. CheiP, A. ChincarinlH A. ChiummJUl H. Ch(B^ 
CN J. Cho\^fil N. ChristenseiP, S. S. Y. ChucP, C. T. Y. Chunpl S. ChunJU, G. CianP, D. E. ClarlJSl J. ClarlEil J. H. ClaytoiCH R Clevid, E. Cocci^P^ R-R CohadoiIH 



- . C. N. ColacinJlS^, J. ColaP, A. Coll^P^ M. ColombinP, A. ContJ^^^ R. ContJH, D. CoolP, T. R. CorbitP, M. CordieiS N. CornislP, A. Cors0, C. A. Cost<P, M. CoughliiP, 
J.-R CouloiP, R CouvareP, D. M. CowarcP, M. Cowar0, D. C. CoynJH, J. D. E. CreightoiP, T. D. CreightoiP, A. M. Cmisjni A. Cumminfl L. Cunninghan0, E. CuocJUl 
R. M. CutleP, K. DahlZEl S. L. DanilishiiP, R. Dannenberfl S. D'AntoniJU, K. Danzmaniltl V. DattilJUl B. Dauder0, H. Daveloz<P, M. DavieJU, E. J. Da\4Sl R. DajP, 
T. DayangJP, R. De RosP^ D. DeBr<P, G. DebreczenP, W. Del Pozz(P, M. del RretJU", T. DenP, V. Dergache^[^ R. DeRos<P, R. DeSalvcR S. Dhurandhai^, L. Di FiorJl, 
^ A. Di Liet&, I. Di FalmflEI M. Di Raolo EmilicP^ A. Di VirgiliJU, M. DiaiH A. Dietfl R DonovaiP, K. L. DoolejP, M. Drag(P^ R. W. R DreveEH J. C. Drigger0, Z. DiP, 
d J.-C. DumaP, T. EberlJZEl M. Edga0, M. EdwardP, A. EffleJ^, R Ehren0, G. EndroczP, R. Enge0, T. Etze0, K. Evan0, M. EvanP, T. Evanfl M. FactouroviclP^I V. Fafon^U^ 
S. FairhursP, Y. FaiP, B. R FarP, D. FazP, H. FehrmanilZil D. FeldbaunP F FeroJ^I. Ferrant&, F Fidecar(P^ L. S. FiniP, I. FiorP R. R FisheP, R. FlaminicP, 
M. FlanigaiP, S. FolejUI E. Fors0, L. A. FortP, N. Fotopoulo0, J.-D. FournieP, J. Fran(P, S. Frasc^P^ F FrasconP, M. FredJzEl M. FreHU Z. FreP, A. FrcisS^ R. FrejIH 
CO T. T. FrickJH, D. FriedriclCEl R FritscheP, V. V. Frolo^il M.-K. FujimotcP, R J. FuldcP, M. FyffJH, J. GaiP, M. GalimbertP, L. GammaitonP^ J. GarciP, F Gam^0^ 

H M. E. GaspaP, G. GemmJH, R. GenJH, E. GeniiP, A. GennaP, L. A. GergeljUl S. GhoslP, J. A. GiaimJlil S. GiampaniP, K. D. Giardin<0, A. Giazott(P, S. GiP, C. Gil0, 
J. GleasoiP, E. GoetflU L. M. GoggiP, G. Gonzale^Il M. L. Gorodetsk>Bll S. GoBleE0 R. GouatjS, C. GraeMl R B. GrafP, M. Granat<P, A. Gran0, S. GraP, C. GrajP, 
I— I N. Gra>0 R. J. S. GreenhalglP, A. M. GretarssoiP, C. GreveriJU, R. Gross(P, H. GrotJzE] S. GmnewalcP, G. M. GuidP^ R. Gupt<P, E. K. Gustafsoi0, R. GustafsoiP, T. H<P, 

OJ. M. HallanP, D. HammeP, G. Hammon(0, J. HankP, C. HanniiP J. Hansoi0, J. HarmP, G. M. HarrjUI I. W. HarrjlSl E. D. HarstacP, M. T. HartmaiP, K. Haughiai0, 
K. HayamJini J.-R HayaiP, J. Heefne0, A. HeidmaniP, M. C. HeintzJUI H. HeitmaniP, R HellcP, M. A. Hendifl I. S. HenJH A. W. Heptonstal0, V. HerrerP, M. HewitsoiCH 
rj S. Hil(0, D. HoalPi K. A. Hodg^, K. Hol0, M. HoltrofP, T. HonJH S. HoopeP, D. J. HoskeiP, J. Hougl0, E. J. HowelP, B. HughejP, S. HusP, S. H. Huttne0, R. IntP, 

T. IsogaP, A. Ivano^[^ K. IzumP, M. Jacobsoi0, E. Jame0, Y. J. Janfi R JaranowskP, E. JessJH W. W. JohnsoiP, D. I. ]onQ0, G. JoneP, R. Jone0, L. jP, R Kalmuffl 
4^ V. KalogerP, S. KandhasamjP, G. KanS J. B. KanneP, R. KasturP, E. KatsavounidiP, W. Katzmai0, H. KaufelMl K. KawabJHl S. KawamuriP, R KawazoflS D. Kelle^P, 
O W. Kell0, D. G. Keppe0, Z. KereszteP, A. KhalaidovskGH F Y. KhalilP, E. A. Khazano^Il! B. KinP, C. KinP, H. KinCH K. KinP, N. KirrP, Y. -M. KinP, R J. Kinjil 

I D. L. Kinze0, J. S. KisseP, S. Klimenk(P, K. Kokeyam<P, V. Kondrasho^[^ S. KorandP, W. Z. Kortt0, 1. KowalskP", D. Kozal0, O. KraniMl V. KringeQ^ S. KrishnamurthjUI 
Q B. KrishnaiP, A. Kr61alP^ G. KuehiQE] R. KumaPl R KweJlE] R K. LanP, M. LandrjP, B. LantP, N. Lastzk<OH C. Lawrifl A. Lazzarin0, R LeacP, C. H. LeJH, H. K. LeJl^ 

H. M. Lejm, J. R. LeoniEl I. LeonoP, N. LerojH*, N. LetendrJlI J. lP, T. G. F L0, N. LiguorP^ R E. Lindquis0, Y. LiiP, Z. LiiP, N. A. LockerbiP, D. Lodhi<P, 
^ M. LorenzinP, V. LoriettP-, M. Lorman(0, G. Losurd(P, J. LouglP, J. LuaiP, M. LubinskP, H. LiiclJSl A. R LundgreiP, E. MacdonalJP, B. MachenschalllZEl M. MacInniP, 
^ D. M. Macleo(P, M. Mageswarai0, K. Mailan(0, E. Majoran<Pf, l. Maksimovi(P, N. MaiP, I. MandeP, V. Mandi(P, M. MantovanP^ A. MarandP, F MarchesonP, 

, ^ , F Marioi0, S. MarkP, Z. MarkJP, A. MarkosyaiP, E. Maro0, J. MarquJUl F MartellP^, I. W. Martii0, R. M. MartiiP, J. N. Mar}0, K. MasoiP, A. Massero0, F Matichar(P, 
L. MatonJU, R. A. MatzneP, N. MavalvalcP, G. MazzolJMl R. McCarthjP, D. E. McClellanP, S. C. McGuirJH, G. McIntyrJH, J. McIveP, D. J. A. McKechaiP, S. McWilliamP, 
G. D. MeadorP, M. MehmelMl T. MeieER A. MelatoP, A. C. MelissinoP, G. MendelP, R. A. MerceP, S. Meshko^[^ C. MessengeP, M. S. Meye0, C. MicheP, L. Milan(0^ 
J. MilleP, Y. MinenkoJlJ, V. R MitrofanoJU, G. MitselmakheP, R. MittlemaiP, O. Miyakaw<P, B. MoP, M. MohaiP, S. D. MohantjP, S. R. R MohapatrP, G. MorencP, 
^> N. Morgad(P, A. MorgiP^ T. MorP, S. R. MorrisP, S. Mosc<0^ K. MossavIZEl B. Mour0, C. M. Mow-Lown|Il C. L. MuelleP, G. MuelleP, S. MukherjeJH, A. MullavejP, 
i/^ H. Miiller-EbhardEH J. MunclP, D. MurphP, R G. Murra}!^, A. MytidiP, T. Nasl0, L. NaticchionP^, V. NeculJP, J. Nelsoi0, G. Newtoi0, T. NguyeP, A. NishizawJP, A. NitJI^ 
F NoceriP, D. NoltinP, M. E. NormandiiP, L. NuttalP, E. OchsneP, J. O'DelP, E. OelkeP, G. H. Ogii0, J. J. OlP, S. H. OlP, B. O'ReilljIH R. O'ShaughnessjP, C. Osthelde0, 
C. D. OtP, D. J. Ottawa^P, R. S. OttenP, H. Overmie0, B. J. OweiP, A. RagiHl G. PagliarolP^ L. PalladincP^ C. Palomb<P, y. PaiP, C. Panko^, F PaolettElHIH 



M. A. PapcEiP M. Paris&^, A. PasqualettP, R. PassaquietP^, D. PassuellcE^, R Pateli( M. Pedrazfl R PeiriP, L. PekowskjP, S. PeniP, A. PerrecJP, G. PersichettP^, 
M. Phelpffl M. RickenpaclJH F Piergiovann&, M. PietkP^, L. PinarJP, I. M. PintcP, M. Pitkii0, H. J. PletsclGH M. V. Pliss0, R. PoggianP^ J. PolJzH F PostiglionjMl 



M. Prat(P, V. PredoP, T. PrestegarcP, l. R. PricJH M. PrijatelPS M. PrincipJ^, S. Priviterfl R. PriJzH G. A. ProdP^ L. G. ProkhoroJH O. PunckeiEEl M. PunturcP, R Pup£cP, 

V. QuetschkJH, R. Quitzow-JameP, F J. RaalP, D. S. RabelinP^ I. RacJU, H. RadkinP, R RaffaP, M. RakhmanoP, B. RankinP, R RapagnanP^ V. RaymoncP, V. R&, 
K. RedwinJH, C. M. ReecP, T. ReecP, T. Re^mbaiP, S. Reic0, D. H. ReitzJUI F RiccP^ R. Riesei0, K. RileP, N. A. RobertsoiHEl F Robine&, C. RobinsoiP, E. L. RobinsoiP, 
A. RocchP, S. RoddJ^ C. RodrigueP, M. RodmclP, L. RollanJP, J. G. Rollin0, J. D. RomancP, R. Roman(0^ J. H. RomiJ^ D. RosinskS^, C. RoveizH S. Rowai0, 
^—H A. RiidigelH R RuggP, K. RyaitnTp SainathaiP, F SalemCEI L. SammuP, V. SandberP V. SannibalJH L. SantamariiH I. Santiago-PrietcR G. SantostasP, B. SassolaP, 

• • B. S. SathyaprakaslP, S. SatcP, p. r. SaulsoiP, R. L. Sava^J^, R- SchillinJZH R. SchnabeuS R. M. S. Schofiel(P, E. SchreibeSH B. SchulJZEl B. F Schutii^ R SchwinberJUl 
^ J. Scot0, S. M. ScotP, F Seifer0, D. Seller0, D. SentenaJUfA. SergeeJH D. A. ShaddoclP, M. ShalteJB B. ShapirJH, R ShawhaiP, D. H. ShoemakeP, A. Siblej^ X. SiemenP, 
• ^ D. Sigfi\ A. Singe0, L. Singe0, A. M. SinteP, G. R. SkeltoiP, B. J. J. SlagmoleP, J. SlutskjP, J. R. SmitP, M. R. Smitl0, R. J. E. SmittP, N. D. Smith-LefebvrJlH K. Somiy^, 
B. Sorazi0, J. Sot(P, R C. Speiritfl L. Sperandi(P^ M. StefszkjP, A. J. SteiiP, L. C. SteiiP, E. SteinerP, J. SteinlechnelZH S. Steinlechneim S. SteplewskP, A. StochinJO, 
^ R. StonJU, K. A. Straii0, S. E. StrigiiP, A. S. StroeeP, R. SturanP^ A. L. Stuve@, T. Z. SummerscaleP, M. SunJU, S. SusmithaP, R J. SuttoiP, B. SwinkelP, M. Tacc<P, 
L. TaffarellcP, D. TalukdeP, D. B. TanneP, S. R TarabriilZEl J. R. TayloIZEl R. Taylo0, R ThomaP, K. A. ThornJCK. S. ThornJH, E. ThranJH, A. Thiirinffl K. V. Tokmako^^H 
C. TomlinsoiP, A. ToncellP^ M. TonellP^ O. Torr&, C. Torre0, C. I. TorriJIH E. Tournefie0, F Travass(P^ G. Traylo0, K. TsenJH, D. UgolinP, H. VahlbmclOEl 
G. Vajent^i^ J. F J. van den Bran(0^ C. Van Den BroeclU, S. van der PutteiP, A. A. van Vegge0, S. Vas0, M. VasutlP, R. VauliiP, M. VavoulidiP, A. VecchicP, G. Vedovat(P, 
J. VeitclP, R J. VeitclP, C. VeltkamfOH D. Verkind0, F Vetran(P^ A. Vicer^2J^ A. E. Villa0, J.-Y. VineP, S. VitalJ^, S. VitalP, H. VoccP, C. VorviclP, S. R VyatchaniiP, 
A. WadP, L. WadP, M. WadiHl s. J. WaldmaiP, L. WallacJH Y. WaiP, M. Wmp\ X. WanJiH, Z. WanJH A. WanneiB R. L. WarP, M. WaP, M. WeinerEH A. J. Weinsteii0, 
R. WeisP, L. WeiEiP R WesselCH M. WesP, T. WestphaOEl K. WettJzEI J. T. WhelaitlTs. E. WhitcomlCP D. J. Whitfi^ B. F WhitinJUI C. WilkinsoiP, R A. Willem0, 

L. WilliamP, R. WilliamJ^R willkJzH L. WinkelmaniQEl W. WinkleOS C. C. WipP, A. G. WisemaiP, H. WittellEl G. Woai0, R. Woole}!^ J. WordeiP, I. Yakushii0, 
H. Yamamot(0, K. YamamotJlil'^, C. C. YancejP, H. YanJ^, D. Yeaton-MasseJH S. YoshidP, R YiP, M. Yver0, A. ZadroznjP, M. ZanoliiP, J.-R ZendrP, F ZhanJH 

L. ZhanJH W. ZhanJH C. ZhacP, N. ZotoP, M. E. ZuckeP, and J. ZweiziJH 
(Affiliations can be found after the references) 

January 16, 2012 



Key words. Gravitational Waves - Methods: observational 



ABSTRACT 



Aims. The detection and measurement of gravitational-waves from coalescing neutron-star binary systems is an important science 
goal for ground-based gravitational- wave detectors. In addition to emitting gravitational- waves at frequencies that span the most 
sensitive bands of the LIGO and Virgo detectors, these sources are also amongst the most likely to produce an electromagnetic coun- 
terpart to the gravitational-wave emission. A joint detection of the gravitational- wave and electromagnetic signals would provide a 
powerful new probe for astronomy. 

Methods. During the period between September 19 and October 20, 2010, the first low-latency search for gravitational- waves from 
binary inspirals in LIGO and Virgo data was conducted. The resulting triggers were sent to electromagnetic observatories for fol- 
lowup. We describe the generation and processing of the low-latency gravitational-wave triggers. The results of the electromagnetic 
image analysis will be described elsewhere. 

Results. Over the course of the science run, three gravitational-wave triggers passed all of the low-latency selection cuts. Of these, 
one was followed up by several of our observational partners. Analysis of the gravitational- wave data leads to an estimated false 
alarm rate of once every 6.4 days, falling far short of the requirement for a detection based solely on gravitational- wave data. 



Key words. Gravitational Waves - Methods: observational 



LSC and Virgo: First Low-Latency LIGO+Virgo Search for Binary Inspirals and Their EM Counterparts 



1. Introduction 

The direct detection and measurement of gravitational-waves 
from coalescing neutron- star and black-hole binaries is a high- 
priority science goal for ground-based gravitational- wave detec- 
tors. Event rate estimates suggest ~ 1/lOOy detectable by the 
initial LIGO-Virgo network rising to ~ 50/y for the advanced 
LIGO-Virgo detector network ( Aba die et al.|201Qa| . Coalescing 
compact binary systems containing at least one neutron star may 
also produce an electromagnetic counterpart to the gravitational- 
wave emission. Observation of both the gravitational and elec- 
tromagnetic emission will allow astronomers to develop a com- 
plete picture of these energetic astronomical events. 

Indeed there is a growing body of evidence that most short, 
hard y-ray bursts (SHGRBs) are the result of binary neutron 
star or neutron star-black hole coalescence (Bloom et al. 2006) 
Berger et al. '^2007, ' Villasenor et"aLl[2Q05t |Berger 2009 ; Nakar 
et al.||2006 ). SHGRBs are known to emit electromagnetic radi- 
ation across the spectrum: Erom prompt, high energy emission 
that lasts ~ minutes, to optical afterglows that decay on the or- 
der of hours and even radio emission that decays on the scale of 
weeks (Nak ar & Pira n"20ir). Coincident electromagnetic and 
gravitational observation of a SHGRB could confirm that the 
central engine of the y-ray burst is indeed a coalescing compact 
binary. 

If the y-ray emission in SHGRBs is beamed, and the cen- 
tral engine is a compact binary merger, there will be many 
compact binary gravitational-wave events for which no y-rays 
are observed by astronomers. On the other hand, optical after- 
glows may be observable further off'-axis than the y-rays al- 
beit somewhat dimmer. There are other proposals for electro- 
magnetic emission from binary neutron- star mergers including 
supernova-like emission due to the radioactive decay of heavy 
elements in the ejecta from the merger ( |Li & Paczy nski 1998 ; 
Metzger et al. 2010 ). This mechanism has been dubbed a "kilo- 
nova" since the luminosity peaks a thousand times higher than a 
typical nova. Kilonova emission is predicted to be isotropic and 
peaks roughly a day after merger. Current blind optical transient 
searches, such as that being undertaken by the Palomar Transient 
Eactory, are only expected to find one of these events per year if 
the neutron binary coalescence rate is at the optimistic end of 
current estimates and that number increases substantially for the 
LSST ( Met zger et al. f2010^. Since the network of gravitational- 
wave detectors is essentially omni-directional, there is a strong 
case to undertake an observing campaign in which compact- 
binary merger candidates, identified in gravitational- wave data, 
are rapidly followed up with electromagnetic observations. 

In this paper we report on the first low-latency search for 
gravitational- waves from compact binary coalescence (CBC) 
in which detection candidates were followed up with optical 
observations. The search covered the period of time between 
September 19 and October 20, 2010, when LIGO detectors at 
Hanford, WA (HI), Livingston, LA (LI) and the Virgo detector 
in Cascina, Italy (VI) were jointly taking data. This time was 
contained within LIGO's sixth science (S6) which began on July 
7, 2009 and Virgo's third science run which began on August 11, 
2010. The group of astronomical partners prepared to make fol- 
lowup observations included the Liverpool telescope, LOEAR, 
the Palomar Transient Eactory, Pi of the Sky, QUEST, ROTSE 
III, Sky Mappe r, Swift, TAROT and the Zadko telescope. A com- 
panion paper ( The LIGO Scientific Collaboration et al. 201 1| ) 
describes more generally the process of collecting gravitational- 
wave triggers from a variety of transient searches, the human 
monitoring of the process and the production of telescope tilings. 



The search resulted in three triggers that passed all of the 
selection cuts. The first occurred during a testing period and no 
alert was sent out. The second trigger was sent out to our astro- 
nomical partners and images were taken by several of them. An 
alert was issued for the third trigger, though its location was too 
close to the sun to be observed. 

This paper is organized as follows: In Sect. II we provide a 
description of the analysis pipeline that produced triggers and 
localized them in the sky. In Sect. Ill we present a performance 
comparison between the low-latency pipeline used to generate 
astronomical alerts and the standard trigger generator used in 
offline searches for compact binary inspirals. Sect. IV presents 
an overview of the run and the details of the triggers that were 
produced. 



2. The Pipeline 

The major components of the analysis pipeline are shown in 
Eig. [T and described in this section. Gravitational- wave data 
from the three participating detectors was first transfered to the 
Virgo site where it was processed by the Multi-Band Template 
Analysis (MBTA) to look for CBC signals. Interesting co- 
incident triggers identified by MBTA were submitted to the 
Gravitational-wave Candidate Event Database (GraCEDb). A 
realtime alert system (LVAlert) was used to communicate infor- 
mation about these events to listening processes that automati- 
cally checked available information for possible anomalous de- 
tector behavior and estimated the location of the source on the 
sky. These automated steps completed within a few minutes after 
data acquisition. 

Eurther processing was needed before an alert was sent to 
our astronomical partners: human monitors reviewed informa- 
tion about the event, consulted the detector control rooms, and 
examined the data quality using a number of higher latency tools. 
Telescope image tilings were generated simultaneously in prepa- 
ration for a positive decision to follow-up a trigger. The entire 
process took 20-40 minutes, with the largest latency incurred 
by the human monitor step. A histogram of the latency incurred 
before the trigger is sent out for further processing is shown in 
Ei g. [2] Details about the rest of the event processing can be found 
in dThe LIGO Scientific Collaboration et al.|20n] ). 



2.1. Trigger production with MBTA 

2.1 .1 . Pipeline structure and parameters 



The Multi-Band Template Analysis (MBTA) ( [Beauville et al. 
2008 ) is a low-latency implementat ion of the standard matched 
filter ( [Wainstein & Zubakov||1962] ) that is commonly used to 
search for gravitational-waves from compact binary coales- 
cences (CBCs). As such, it relies on the accurate mathematical 
models of the expected signals to be used as search templates. 
Time domain templates with phase evolution accurate to sec- 
ond post-Newtonian order were used in the search. The search 
covered sources with component mass between 1-34 M© and 
total mass below 35 M©. A fixed bank of templates constructed 
with a minimal match of 97% was used to scan this parame- 
ter space. The template bank was constructed using a reference 
noise power spectrum taken at a time when the detectors were 
performing well. Given the detectors' typical noise spectra, the 
low-frequency cutoff' of the templates was set to 50 Hz, thus 
keeping the computational cost of the search light while losing 
negligible signal-to-noise ratio (SNR). 
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match within expected uncertainties. The mass coincidence cri- 
terion is 



Sky Localization 
and 

Data Quality Check 



Fig. 1. An overview of the pipeline. Data produced at each of the 
three detector sites is transfered to a computer at Cascina where 
triggers are produced by MBTA and sent to GraCEDb for stor- 
age. Upon receiving events, GraCEDb alerts the sky localization 
and data quality check processes to begin and the results are then 
sent back to GraCEDb. If the triggers are localizable and of ac- 
ceptable data quality then they are sent out for further processing 
and possibly followup by an optical telescope. The double stroke 
connections in the diagram are provided by LVAlert. 



An original feature of MBTA is that it divides the matched 
filter across two frequency bands. This results in two immedi- 
ate benefits: (1) The phase of the signal is tracked over fewer 
cycles meaning sparser template banks are needed in each fre- 
quency band and (2) a reduced sampling rate can be used for the 
lower frequency band, reducing the computational cost of the 
fast Fourier transforms involved in the filtering. The full band 
signal-to-noise ratio p is computed by coherently combining the 
matched filtering outputs from the two frequency bands. The 
boundary between the low and high frequency bands is chosen 
such that the SNR is roughly equally shared between the two 
bands. The boundary frequencies range from 110 Hz to 156 Hz, 
depending on the detector and on the mass ratio of the template 
bank of the individual MBTA processes. 

A trigger is registered when p > 5.5. Triggers are clus- 
tered across the template bank and across time: triggers less than 
20ms apart are recursively clustered under the loudest trigger. 
Clustered triggers are subjected to slx^ test to check if the SNR 
distribution across the two frequency bands is consistent with the 
expected signal. The discriminating power of such a 2-band 
test is not as high as in typical implementations based on a larger 
number of frequency bands, but off'ers the advantage of having a 
negligible computational cost in the multi-band framework. The 
test can therefore be applied to all triggers at the single detector 
level. Triggers pass the;^^-test if 



;^2< 3 (2 + 0.025p2). 



(1) 



Single detector triggers that pass the test are tested for 
coincidence across detectors. Triggers from two detectors / and 
j are considered coincident if their time and mass parameters 



\Mi-Mj\ 



0.05 



Mi + M, 



(2) 



where the chirp mass is given by At = (mim2)^^^(mi -\- m2)~^^^ 
in terms of the component masses, mi and m2. The time coinci- 
dence criterion is 



1^/ -tj\ < Atij 



(3) 



where the maximum allowed time delays AtniLi = 20 ms and 
Atmvi - AtLiYi - 40 ms account for both the time of flight 
of the signal from one detector to another and for the exper- 
imental uncertainty in the arrival time measurement. The ar- 
rival time is measured by performing a quadratic fit around the 
maximum of the match filtered signal. It is then extrapolated 
to the time when the gravitational-wave signal is at a refer- 
ence frequency chosen to minimize the statistical uncertainty 
on the measurement. This has the important eff'ect of reducing 
the background by allowing tighter coincidence and improv- 
ing the accuracy of positio n reconstruction by tr iangulation. It 
has been shown elsewhere ( Acemese et al.||2007| ) that the opti- 
mal reference frequency is approximately that frequency where 
the detectors are most sensitive, although it depends somewhat 
on the mass of the source. A detailed study of simulated sig- 
nals in S6A^SR3 noise resulted in the empirical parametrization 

/reference = [170 - M(5.1/Mo)] Hz. 

Triple detector coincidences are identified as a pair of HlLl 
and HI VI coincidences sharing the same HI trigger. Although 
the pipeline identifies both double and triple coincidences, only 
the latter were submitted to GraCEDb, in line with the intention 
to focus on remarkable events, for which a sky location could 
be extracted by simple triangulation. Each triple has associated 
with it a combined SNR pc defined by 



7G{H1,L1,V1] 



(4) 



and a false alarm rate (FAR). Triggers with a smaller FAR, i.e. 
a larger p^, are more likely to be gravitational- waves. The non- 
stationary nature of the background means that a simple map- 
ping from a given pc to a corresponding FAR does not exist. 
Instead, the FAR must be explicitly calculated for each trigger. 

The FAR is estimated as follows. Let Ti be the analysis time 
needed to accumulate the last 100 triggers in detector /. Let 
^louder be the number of mass-coincidences that can be formed 
from the last 100 triggers in each detector for which pc is greater 
than the combined SNR of the observed coincidence. Then the 
expected rate of coincident triggers arising from background 
may be estimated as the product of the individual detector rates 
multiplied by a factor that accounts for the coincidence windows 
in time and mass. In particular, we define the FAR associated 
with a triple-coincident trigger to be 

A Atmii AtHivi 

¥AR = a Mouder^^^ 7^^^ 

^ HI ^ LI ^Vl 



(5) 



where a is an empirical factor with a value of 2 tuned to adjust 
the estimated FAR to the average observed rate of such triple- 
coincident triggers. 

The low latency search is dependent on the short-order avail- 
ability of the strain time series h(t) measured by the detectors. 
Tools to reconstruct h(t) with very low latency have been devel- 
oped over the last few years, and had reached a mature state by 
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the time of the S6/VSR3 runs. Calibration accuracies better than 
15% on the amplitude of h(t) were typically achieved for all de- 
tectors, which is quite sufficient for the purposes of our search. 

For each detector, the h(t) channel is produced at the ex- 
perimental site computing center. The HI and LI h(t) data are 
then transfered to the Virgo site computing center, where all the 
MBTA processing takes place. To minimize the latency, the com- 
munication of input/output data between the different processes 
involves no files on disk. It relies instead on the use of shared 
memories and of a TCP-IP based communication protocol de- 
veloped for the Virgo data acquisition system. The whole set of 
processes ran on six computers. 

2.2. GraCEDb and LVAIert 

The storage and communication capabilities of the pipeline in 
Fig. [T] are provided by the Gravitational- wave Candidate Event 
Database (GraCEDb) and the LValert messaging system. The 
purpose of these technologies is to provide a continuously run- 
ning system to ingest, archive and respond to gravitational-wave 
triggers. Communication with individual telescopes is handled 
at a later stage and uses whatever protocol is appropriate for the 
particular telescope. 

The GraCEDb service stands b ehind 
[Apache Software Foundation] |2Q 11 1) 



an Apache 
server and is built 



(The 



on 



Django (Django Software Foundation! |2Q11|, a command line 



client that uses an HTTP/ReST ( [Fieldingl |2000l [Fielding 



|et al.[[2QQ2]) interface to the server and authenticates with 



X5Q9 ([Housley et al.[ri999[ ) certificates to a MySQL C Oracle 
[Foundation[[2010[ ) database back-end. A command-line client 
allows easy automation of event submission to GraCEDb. The 
prototype system used during S6A^SR3 was capable of ingest- 
ing triggers from MBTA and a number of other search pipelines. 
The raw trigger information was stored in an easily accessible 
archive and the most relevant information about the trigger, such 
as the time and significance, were ingested into the database. 
Upon successful ingestion, the trigger is given a unique identi- 
fier that is returned to the submitter. An alert is then sent out via 
LVAIert. 

LVAIert is a communication client built on XMPP dThel 



XMPP Standards Foundation 2011) technology and the PubSub 



extension (Millard et al. 201 1 ). The system allows users to create 
nodes to which information can be published; users subscribe to 
nodes from which they want to receive that information. In the 
context of the current search, a node was set up for sending out 
alerts about MBTA triggers. A listener client is also provided 
that waits for alerts from the nodes to which the user is sub- 
scribed. By default, the listener simply prints any information it 
receives to the stdout, but it also allows users to develop plu- 
gins which can take action in response to the alerts. This is the 
mechanism for launching the data quality and sky localization 
jobs in Fig.[T] 

2.3. Event processing 

While MBTA is responsible for producing coincident triggers, 
there remains further processing to determine whether or not 
the triggers are suitable for external followups. In particular, this 
processing consists of performing sky localization, checking the 
quality of the data and determining whether or not the event is a 
known hardware injection, as described below. The mechanism 
for performing these tasks is an LVAIert listener that responds to 
alerts sent out by GraCEDb in response to new MBTA triggers. 



Time after the event (minutes) 

Fig. 2. Total latency of the automatic processes during the 
S6/VSR3 run. The x-axis is the diff'erence between the time at 
which the trigger was available for human monitoring and the 
trigger's time. This includes an average of 63 s for all the data to 
be available at the Virgo site computing center and an average 
of 142 s for the trigger to be submitted to GraCEDB. This does 
not include the time for the human monitoring (~20-40 minutes) 
that takes place before an event is sent out for electromagnetic 
followup. 

The sky localization procedure proceeds by computing, as a 
function of location on the sky for a pre-determined grid, the 
errors in timing and amplitude. Probabilities are assigned by 
comparing these measured quantities with distributions of these 
quantities obtained from simulated signals. More specifically, 
the measure of timing accuracy used is 



ti)-(tj(a,S)-tj)r 



(6) 



where ti(a,6) is the geocentric arrival tim^of a source located 
at right ascension a and declination S in detector / and ?/ is the 
measured arrival time. To improve sky localization performance, 
the time of arrival is taken to be the time the signal crosses the 
reference frequency ( Acernese et al.[[2007, ) of 140Hz. The am- 
plitude consistency parameter is given by 



/r)(02 



(7) 



where the (/) superscripts label the detector, Dq^ is the effective 
distance and we have suppressed the (a, S) of every quantity on 
the right hand side of the equation. Analytically, the eff'ective 
distance is defined by 



1 + COS^ L 



+ fI cos^ l 



-1/2 



(8) 



where D is the actual distance to the source, x = F+^x(o:, <A) 
are the response functions of the detector (that depend on the po- 
larization angle, i/^) and l is the inclination angle of the source, 
and the (/) superscript is understood. In practice, the matched 
filtering procedure only provides a measurement Dq^ that is not 
typically invertible to obtain Z), F+^x, etc. We introduce the quan- 
tity (suppressing the (/) superscript once again) 



1 



(9) 



We assume that gravitational- waves travel at the speed of light. 
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to provide an ad hoc measure for determining the consistency 
of the amplitudes measured in each detector with a particular 
sky location. In practice, the use of this quantity improves the 
sky localization accuracy by roughly a factor of two, mean- 
ing it helps to break the mirror-image degeneracy inherent in 
a three-detector network when only triangulation is used to lo- 
cate the source on the sky. Values of A^rss(^^, ^) and AQyss(o;,S) 
are assigned probabilities by comparing their values to those 
of a distribution of simulated signals. More specifically, simu- 
lated signals are placed into detector data and probability dis- 
tributions are computed from the values of A^rss(Q^tme, ^tme) and 
^Qrssio^tmQ^ ^tme). where the subscript indicates that the quantity 
is evaluated at the true location of the source in the sky. Since 
there is negligible correlation between the timing and amplitude 
measures, the normalized product of their values yields the prob- 
ability distribution over the sky, i.e., the skymap. 

As a demonstration of the performance of the sky localiza- 
tion routine, 10122 simulated signals were injected into data 
taken from the 6^^ week of S6/VSR2. The locations on the 



sky and distances were given by a galaxy catalog ( Kopparapu 



|et al. 1200 8]). More details about these injections can be found in 
Sect. |3.2[ The SNR distribution of the population of injections 
recovered by the pipeline is given in Fig. |3] We use two figures 
of merit to assess sky localization performance. The first is a 
quantity we call searched area. It is the area on the sky con- 
tained in the pixels that are assigned a probability greater than 
the probability assigned to the pixel containing the true source 
location. In other words, it is the smallest area on the sky that 
one would be required to image to ensure the true source loca- 
tion is imaged. A complimentary figure of merit is the angular 
distance between the true source location and the most probable 
location identified by the sky localization routine. These figures 
of merit are shown as cumulative histograms in Fig. |4] To im- 
prove the performance of the sky localization routine, a galaxy 
catalog is used to guide the pointing in the following way: The 
skymap is weighted by the fraction of the cumulative blue lumi- 
nosity (used in this context as a proxy for star formation rates) in 
each pixel out to the smallest eff'ective distance measured by any 
of the detectors ( The LIGO Sc ientific Collaboration et al.|201 Ij ). 
The eff'ect of using a galaxy catalog is evident in the figures of 
merit shown in the upper curves in Fig.|4] Focusing on the lower 
curves, the angular distance figure of merit shows that for 60% 
of the injections the maximum likelihood point on the skymap 
is < 5 degrees away from the true source location, whereas the 
searched area is a much larger region on the sky. This indicates 
that although the maxima of the skymaps are close to the true 
source locations, they are quite broad. Application of the galaxy 
catalog does not do much to the angular distance histogram, 
but it eliminates large regions of the sky, eff'ectively making the 
broad peaks much sharper. In Fig. |5] the median searched area 
is plotted as a function of combined SNR. The upper panel de- 
picts the results without the aid of the galaxy catalog, while the 
lower panel includes those improvements. A least squares fit of 
the median searched areas is given by the red lines. The fit is 
to a functional form ao + <2_i/p + a-^lp^ . The results that make 
use of the galaxy catalog should be understood as upper limits 
on th e im provements a galaxy catalog can bring. As described in 
Sect. 3^ the injection distribution was such that injections were 
more likely to come from galaxies with larger blue light lumi- 
nosity. Although the choice to use the blue light luminosity is 
well motivated, it does not include, for example, binaries arising 
in elliptical galaxies, or the possibility that binaries follow some 
other property of their host galaxies. Better astrophysical priors 
are a subject we leave for future investigations. 




Combined SNR 



Fig. 3. The combined SNR distribution of the family of found 
injections used to demonstrate the sky localization performance. 
Emphasis was placed on weaker signals, where a first detection 
is more likely to arise. 
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Fig. 4. Sky localization performance with and without the use of 
a galaxy catalog. The upper pane shows a cumulative histogram 
of the searched area in square degrees. The lower pane is a cu- 
mulative histogram of the angular distance, in degrees, between 
the injected location and the maximum likelihood recovered lo- 
cation. In both plots the red solid line is the performance with 
the aid of the galaxy catalog and the blue dotted line is the per- 
formance without the galaxy catalog. 



P combined 




Fig. 5. Sky localization performance as a function of combined 
signal-to-noise ratio with and without the use of a galaxy catalog 
(bottom and top panel, respectively). The blue bars indicated the 
median searched area, in square degrees, in each bin and the red 
lines depict a fit to these values. 
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Data quality plays an important role in offline searches, 
where detailed studies are performed to select data free from 
obvious instrumental or environmental artifacts by using infor- 
mation in auxiliary channels. While all the data quality infor- 
mation is not available with low latency, online production of 
a significant number of data quality flags was implemented dur- 
ing S6/VSR3. This allows for quick feedback on the current data 
quality, and to check the data at the time of a trigger against a list 
of trusted data quality flags produced online. As new data quality 
flags are developed and tested, the list is updated accordingly. 

If an event passes all of our checks then it is sent out to 
LUMIN, where it is then processed further and possibly sent out 
for electromagnetic observation. 



3. MBTA Performance and Validation 

In this section we look at the performance of MBTA. There 
are typically two ways in which performance is assessed: (1) 
through hardware injections, where the mirrors of each detec- 
tor are physically displaced to simulate the presence of a signal 
and (2) software injections in which simulated signals are placed 
into detector data. We will first look at MBTA's performance on 
hardware injections and then provide a detailed co mparison of 
its performance relative to the CBC offline pipeline ( A bbott et al. 
[20081 [20071 [20Q9a|bt[AFadie al.|201Qb| ), known as iHope, on 
a set of common software injections. 



o 



3.1. Hardware injections 

Hardware injections are produced through the displacement of 
the mirror located at the end of one of the detector's arm. They 
are performed coherently in the three detectors, i.e. with injec- 
tion time, amplitude and phase chosen in each detector to be 
consistent with the location of the source on the sky. During the 
S6/VSR3 scientific run, there were approximatively three hard- 
ware injections per day in each detector simulating the coales- 
cence of a compact binary system, with a period of intensive 
injections between August 27, 2010 and September 3, 2010. For 
a number of reasons, not every injection was successful at every 
detector. 

Two families of non-spinning waveforms were used for 
the hardware injections. The first family corresponds to 
analytic inspiral-merger-ringdown waveforms based on the 
Eff'ective One-Body (BOB) model extended and tuned to match 
Numerical Relativity simulations of binary black hole coales- 
cences i Buonanno et al.|2007][Pan et al.|2008| ). The second fam- 
ily of waveforms corresponds to restricted parameterized 2PN 
inspiral waveform computed in the time domain (Blanchet 19961 
[Arun et SI1[2004| ). Among the hardware injections successfully 
performed in all three detectors, 62% are from the first family 
and 38% are from the second. 

The parameters of the hardware injections were adjusted to 
appear in each detectors with an SNR < 100. The total masses 
of the simulated binary systems were distributed between 2.8 M© 
and 35 M©, with the masses of the components between 1.4 M© 
and 35 M©. The hardware injections were distributed between 
1 Mpc and 80 Mpc, and the sky locations were chosen randomly. 
Figure [6] shows the chirp mass distribution of all the hardware 
injections which were successfully found by MBTA during the 
S6/VSR3 run. 

An injection is considered to be recovered if MBTA pro- 
duced a trigger with an end time within 20 ms of the end time 
of the injection. We restrict our attention to injections occurring 



In 



Mchirp(M^ 



Fig. 6. The chirp mass distribution of all the hardware injections 
found in a single detector by both MBTA and the offline iHope 
pipeline during the S6/VSR3 run. 



when all three detectors are taking science quality data and all 
of the MBTA processes are up and running. Our focus is further 
narrowed by the requirement that the data quality is acceptable in 
60 s preceding a trigger. This is to ensure that the Fourier trans- 
form required by the matched filter can be performed. Given 
these caveats, MBTA achieved a 93% efflciency in detecting 
hardware injections in triple coincidence during the S6/VSR3 
run. Figure [7] shows these recovered hardware injections and in- 
dicates the type of coincidence (double or triple) they were found 
in. 
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Fig. 7. This plot shows all the successful hardware injection 
which took place coherently in the three detectors. 93% of those 
which where above the threshold (expected SNR > 5.5) in each 
detector were detected in triple coincidence by MBTA. 



An important step in the validation of MBTA is the com- 
parison of the injected signal parameters with those recovered 
by MBTA. Figure [8] shows that the parameters of the injected 
signals are in good agreement with those recovered by MBTA. 
The first plot of Fig. [8] shows the timing accuracy parameter of 
the LIGO- Virgo network, A^rss, given by Eq.[6] The second and 
third plots of Fig. [8] show the diff'erence between the chirp mass 
and eff'ective distance of the binary injected and observed. From 
these three plots, it is clear that the low-latency search provides 
parameter estimations with comparable accuracy to that of the 
offline iHope pipeline. 
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Fig. 8. These plots show the differences between the parameters 
of the hardware injections and those of the triggers found by 
both MBTA and the offline iHope pipeline. The first plot shows 
the timing accuracy parameter A^rss of all triple coincidences. 
The second plot corresponds to the diff'erence between the chirp 
mass of the binary injected and observed. The third one corre- 
spond to the diff'erence between the eff'ective distance of the bi- 
nary injected and reconstructed by MBTA. The second and third 
plots contain all the hardware injections observed in HI, LI or 
VL 



3.2. Software Injections 

Hardware injections are for obvious reasons impractical for large 
scale performance studies. Instead, it is preferable to inject simu- 
lated signals into data that is already on disk. Here we will com- 
pare the performance of MBTA with the existing offline CBC 
pipeline, known as iHope. Both pipelines were run on a week of 
data containing simulated binary neutron star and neutron star- 
black hole binary signals. For this study a total of 10122 soft- 
ware injections were performed and analyzed in multiple paral- 
lel runs. 



Simulated waveforms we re generated using the stationary 
phase approximation (SPA) (Droz et al. 1999[ Thome 1987[ 
|Sathyaprakash & Dhurandhar||1991| ) with the amplitude ex- 
panded to Newtonian order and the phase expanded to sec- 
ond post-Newtonian order with the upper cut-off' frequency 
at the Schwarzschild innermost stable circular orbit (ISCO). 
Their locations were randomly chosen from the Compact Binary 
Coalescence Galaxy Catalogue ( Kopparapu et al.|2 QQ8 ) in such 
a manner that the probability of choosing a galaxy is propor- 
tional to its blue light luminosity. Injections were made out to 
a distance of 40Mpc which roughly corresponds to the limits 
of sensitivity of the LIGO-Virgo network during the observa- 
tion period, for the considered sources. The mass of each bi- 
nary component ranged from 1 M© to 15 M© with the constraint 
on the total mass of the binary to be between 2 M© and 20 M©. 
These choices were made to focus on systems likely to contain 
at least one neutron star, where an electromagnetic counterpart 
is expected. 

After data quality vetoes were applied MBTA and iHope re- 
covered diff'erent numbers of triple-coincident triggers, 736 and 
859, respectively. Both pipelines recovered 709 identical injec- 
tions in triple coincidence. Among the rest of the MBTA triples, 
twenty five were found by the offline CBC pipeline in double 
coincidence and the remaining two were completely missed. All 
these signals were near the threshold of detectability in one or 
more detectors. 



Timing accuracy The primary goal of the low-latency pipeline 
is to send triggers out to the astronomical community for elec- 
tromagnetic followup observations. Hence, localizing the GW 
candidate event on the sky is one of the essential parts of the 
search. Good timing accuracy for recovered injections is essen- 
tial for good sky localization. Figure [9] shows the normalized 
distributions of the timing accuracy parameter of the LIGO- 
Virgo network, A^rss, given by (|6]). From this plot it is clear that, 
overall, MBTA's performance in recovering the arrival time of 
gravitational- wave is better or comparable to that of the offline 
search. This is expected because MBTA uses a quadratic fit to 
find the peak of the SNR time series whereas the offline pipeline 
simply takes the maximum of the time series. 
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Fig. 9. Normalized distribution of timing accuracy of triggers 
detected by MBTA and CBC offline analysis pipeline. MBTA 
shows slightly better performance than the offline analysis. 



Efficiency as a function of distance is another key characteris- 
tic of an analysis pipeline. For each pipeline we measure its ef- 
ficiency at recovering a GW signal as a triply coincident trigger. 
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Fig. 10. Detection efficiency of MBTA and CBC offiine analysis 
pipeline as a function of distance. Shaded regions indicate un- 
certainties. MBTA detects all simulated gravitational- wave sig- 
nals injected in the nearby universe that the offiine analysis does, 
but finds systematically less number of signals than the offiine 
pipeline beyond D > 3 Mpc. 



requiring also a false alarm rate less than or equal to 0.25 events 



per day, as for the alert generation. The result is shown in Fig. [TO 
At distances of about 3 Mpc, MBTA begins to systematically re- 
cover fewer signals than iHope. As expected this trend continues 
as both pipelines lose efficiency at large distances. Much of the 
decrease in efficiency of MBTA can be attributed to the fact that 
it imposes an eff'ective threshold at a slightly higher SNR for a 
signal to be detected. 

Chirp mass and SNR. Finally, Fig. [TT] shows the SNRs and 
chirp masses in HI recovered by the pipelines. They are in very 
good agreement with each other. The sparse density of the tem- 
plate bank at high mass shows up as some discreteness in the 
values of chirp mass recovered by MBTA, because MBTA uses 
a constant template bank, whereas iHope recomputes the noise 
power spectrum every 2048 s and a new template bank is pro- 
duced in response. Similar agreement is found for the other two 
detectors. 

In summary, the performance of MBTA is comparable to that 
of the offiine pipeline. 



4. Results 

In this section we present the results of the analysis. We begin 
with an overview of the joint LIGO- Virgo science run before 
turning our attention to the triggers that passed the selection cuts. 



4.1. The S6/VSR3 run 

The joint LIGO/Virgo data taking started on August 11, 2010 
when Virgo started its 3rd science run (VSR3), joining LIGO's 
6th science run (S6). After a test and adjustment period during 
the first part of the S6/VSR3 run, the software of the low latency 
pipeline was frozen on August 27. From this time on, until the 
end of the run on October 19, the full pipeline was operating 
in production mode. The trigger production was first monitored 
to validate the pipeline during online operations and then the 
submission of EM alert was enabled on September 19. 

During the production period of 52.6 days, the science mode 
duty cycle of the HI, LI and VI detectors were respectively 
63.9%, 64.8% and 69.7%. The three detectors were operating 
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Fig. 11. Comparison of chirp mass and SNR between MBTA and 
CBC offiine analysis pipeline. 



simultaneously in science mode for a total time of 18.2 days, 
corresponding to a 34.6% duty cycle. 

The duty cycle of MBTA during the triple coincidence time 
was 94.2%, 98.7% and 97.8% for the single HI, LI and VI trig- 
gers and 91.2% generating triggers in triple coincidence. Most 
of the down time occurred in one of a small number of periods 
lasting a few hours. The main problem was temporary network 
overloads between the LIGO Hanford site and Virgo Cascina 
site. The resulting delay in the arrival of HI data prevented it 
from being used by MBTA. 

Over the S6/VSR3 production period, 89 triple coinci- 
dences - including hardware injections - were detected by the 
MBTA pipeline and submitted to the GraCEDb database. A 
few of them (10) triggered multiple submissions of the same 
loud event corresponding to nearby "satellite" events. Three 
GraCEDb submissions failed. One because of a GraCEDb disk 
access problem. The other two failures were due to a problem 
with network authentication. 

Future operations with improved configurations, software 
versions and monitoring tools are expected to reduce the down 
time of the pipeline which involved five diff'erent computing fa- 
cilities located in Europe and North America. 
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4.2. Triggers 

After removing the hardware injections, a total of 42 triple coin- 
cident triggers were observed during the search. The application 
of the online data quality flags reduced the number of triple coin- 
cident triggers from 42 to 37. The time coincidence window was 
chosen conservatively (larger than the light travel time between 
sites) and, as a result, only 23 of these triggers were localizable 
on the sky. 

At this stage of the pipeline, the triggers were passed to 
LUMIN which generated alerts for events having a false alarm 
rate of less than 0.25 events per day (1 event per day up to August 
31). This cut reduced the number of triggers to 13 which were 
passed to the control rooms and on-call experts for further qual- 
ity assessments. 

Out of these 13 possible alerts generated by LUMIN, only 
3 met the requirement of having at least one neutron star (M < 
3.5 Mq) associated with the merger. Table [T] gives a snapshot of 
the parameters of these three triggers. 

The first trigger, G 16901, from August 31, occurred during 
an initial testing period before alerts were sent out. A decision 
on this trigger was reached in the control room 14 minutes af- 
ter the trigger time. The second trigger, G20190, on September 
19, was accepted and images of the corresponding sky location 
were taken by Quest, ROTSE, SkyMapper, TAROT and Zadko. 
The decision to issue an alert was reached 39 minutes after the 
event occurred. The image analysis is in progress. Eigure 12 
shows the skymap produced for this trigger. The 90% confi- 
dence region was reduced from nearly 600 square degrees to 
3.3 square degrees with the application of the galaxy catalog. 
The third trigger, G23201, on October 6, was unfortunately lo- 
cated too close to the sun, making it impossible to image. The 
decision to send the trigger out occurred 16 minutes after the 
trigger. Overall, these three triggers have SNR values close to 
the threshold value, with a false alarm rate of one per few days, 
typical of the expected background triggers. They are therefore 
not detection candidates on the basis of the gravitational-wave 
data analysis. 




0.0000 0.0015 0.0030 0.0045 0.0060 
Probability per deg^ 



0.0075 0.0090 



Fig. 12. Skymap for the G20190 trigger on September 19. 



Eigure 13 shows the cumulative rate of observed triggers as 
a function of the upper threshold applied on the estimated false 
alarm rate. The distribution focuses on triggers collected during 
the production period with a EAR less than 200 per year, requir- 
ing at least one neutron star (M < 3.5 M©) and excluding hard- 
ware injections. The vertical line indicates the threshold of EAR 



lower than -91 per year (0.25 per day), which was used to de- 
termine which triggers were candidates for EM follow-up. This 
figure shows that the online EAR estimation is reasonable and 
therefore the background is under control. In particular, this fig- 
ure shows no evidence of the EAR being underestimated, which 
is important since we do not want to unduly promote uninterest- 
ing triggers. 



Uppeit threshold; 



03 
> 
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False alarm rate (1/yr) 

Fig. 13. Cumulative rate of triggers observed during the produc- 
tion period (excluding hardware injections) as a function of the 
upper threshold applied on the estimated false alarm rate. 



5. Discussion 

The coalescence of binary systems containing neutron stars is 
the most promising source for the detection of both gravitational 
and electromagnetic radiation. We have presented the first low- 
latency search for the gravitational-waves from these systems 
during the period between September 19 and October 20, 2010. 
The search resulted in a single trigger with a false alarm rate of 
one per 6.4 days being followed up with optical telescopes. The 
results of the image analysis are pending. 

This exercise has resulted in a low-latency search for binary 
inspirals that performs on par with the standard offline pipeline. 
Triggers are produced on the scale of minutes and it is likely 
this can be further reduced. The limiting latency in sending trig- 
gers to electromagnetic observatories is the human monitoring 
involved. However, our demonstration that reliable false alarm 
rates can be computed rapidly suggests that this step could be 
removed in the advanced detector era. 

Improvements of the pipeline will be explored in the future. 
Eor instance, the volume probed could be optimized by apply- 
ing thresholds which could depend on the sensitivity and on the 
data quality of each detector. More detailed data quality infor- 
mation could help for this step. The search area in the sky might 
be reduced by using some coherent technique, possibly in a hier- 
archical way. Extending the emission of fast alerts to significant 
double coincidences may also be useful. 

Advanced LIGO and Virgo are likely to detect -50 neutron 
star coalescences per year ( Abadie et al.|2010a| . Successful ob- 
servation of joint EM-hGW emission depends crucially on us- 



10 



LSC and Virgo: First Low-Latency LIGO+Virgo Search for Binary Inspirals and Their EM Counterparts 
Table 1. Parameters of the three triggers which passed all the selection cuts. See text for details. 



Detector SNR D^s[Mpc] mdMe] m2[Me] M[Mq] 



G169Q1: 967254112; Combined SNR=9.99; FAR-^^1.1 days 



HI 6.15 55 1.03 
LI 5.61 54 1.36 
VI 5.52 19 1.35 


2.06 
1.38 
1.37 


1.26 
1.19 
1.18 


G20190: 968932960; Combined SNR= 


10.0; FAR- 


"^=6.4 days 


HI 6.07 99 2.94 
LI 5.65 106 3.05 
VI 5.60 27 2.23 


3.00 
3.11 
4.15 


2.59 
2.68 
2.62 


G23201: 970399241; Combined SNR= 


10.1; FAR- 


"^=5.5 days 


HI 5.75 58 1.04 
LI 5.84 41 0.98 
VI 5.96 14 0.97 


1.99 
1.95 
1.91 


1.24 
1.19 
1.17 



ing all the available information on these sources. Expectations 
for electromagnetic emission and their lightcurves, for example, 
will be important for designing optimal observing campaigns. 
In addition, even with a three-detector network, the sky local- 
ization ability is limited. More detectors, such as the LCGT or 
some other future detector will be important for increasing the 
chances of making successful joint EM-fGW observations. 
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